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Abstract. The implementation of fault-tolerant quantum gates on encoded 
logic qubits is considered. It is shown that transversal implementation of 
logic gates based on simple geometric control ideas is problematic for realistic 
physical systems suffering from imperfections such as qubit inhomogeneity or 
uncontrollable interactions between qubits. However, this problem can be 
overcome by formulating the task as an optimal control problem and designing 
efficient algorithms to solve it. In particular, we can find solutions that implement 
all of the elementary logic gates in a fixed amount of time with limited control 
resources for the five-qubit stabilizer code. Most importantly, logic gates that 
are extremely difficult to implement using conventional techniques even for ideal 
systems, such as the T-gate for the five-qubit stabilizer code, do not appear to 
pose a problem for optimal control. 
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1. Introduction 

Quantum information processing has been a topic of intense theoretical and 
experimental research for several years. Many potential physical realizations of 
quantum computers have been proposed, and although the experimental realization 
of quantum information processing remains a challenge, there have been many 
experimental accomplishments including the demonstration of control of multi-qubit 
dynamics in spin systems using nuclear magnetic resonance techniques [1] and all- 
optical quantum information processing using photons [2J, for instance. While liquid- 
state NMR and optical quantum computing may suffer from inherent scalability issues, 
there has also been significant progress in ion-trap architectures and several proposals 
for making such architectures scalable exist 011]. In solid-state systems successes have 
been more modest but controlled interactions of quantum dots have been demonstrated 
in some systems [5j |6] . 

One major obstacle to scalability is the increasing difficulty in effectively 
controlling the dynamics of many qubits as the complexity of experimental systems 
increases. While it is in principle easy to control a single two-level system, the 
situation is more complicated when many qubits are involved. In many cases individual 
addressability of qubits in a register or array is difficult to achieve; inhomogcneity 
may result in different qubits having different responses to control fields, and the 
dynamics is complicated by uncontrollable inter-qubit couplings. All of these issues 
present challenges, even in the absence of environmental noise and decoherence. These 
problems are magnified because the necessity to be able to correct at least a certain 
amount of inevitable errors means that bits of quantum information must be encoded 
in logical qubits consisting of multiple physical qubits to ensure that there is sufficient 
redundancy. Even if only the most basic level of protection is assumed, at least five 
physical qubits are required to encode a single logical qubit so that we can correct a 
single (bit or phase flip) error |10j . In practice multiple layers of error correction will be 
necessary to achieve fault-tolerant operation of a quantum processor for a sufficiently 
long time to allow the completion of a non-trivial quantum algorithm |10j . Thus, in the 
setting of fault-tolerant computation using encoded qubits even the implementation 
of single qubit logic gates becomes a non-trivial multi-qubit control problem. 

The implementation of fault-tolerant gates quantum logic gates seems to be a 
problem that is well suited for optimal control. Optimal control theory has been used 
in various papers to implement two and three qubit gates as well as simple quantum 
circuits, and has been shown to improve gate fidelities, gate operation times and 
robustness (e.g. [7HHJG3). ^ n this paper we consider the implementation of logic gates 
on encoded qubits, which is a next logical step. Although most logic gates on encoded 
qubits for the most common codes can be implemented, in principle, by applying a 
particular single qubit operation to all physical qubits, what is known as transversal 
implementation, not all logic gates can be implemented this way even in the ideal case, 
and transversal implementation is problematic in the presence of inhomogcneity or 
uncontrollable couplings between physical qubits, and the implementation of encoded 
logic gates in this setting is therefore in an interesting challenge for optimal control. 
In particular we study the implementation of a set of logic gates for the five-qubit 
stabilizer code for model systems with both inhomogeneity and fixed inter-qubit 
couplings. The paper is organized as follows. In Sec. [2] we briefly introduce stabilizer 
codes and the construction of fault-tolerant logic gates with emphasis on the five-qubit 
stabilizer code that will be the main focus of this paper. In Sec. [3] we describe two 
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types of model systems and discuss the implementation of encoded logic gates for these 
system based on geometric control ideas. In section|4]we formulate the optimal control 
problem and show how to design iterative algorithms to solve it numerically. In Sec. [5] 
the results are applied to find optimal controls to implement a set of standard logic 
gates for the five-qubit stabilizer code for two classes of model systems, and the merits 
and drawback of the optimal control approach compared with the simpler geometric 
control schemes are discussed. 

2. Stabilizer codes and fault-tolerant gates 

The state of a single qubit can be represented by a density operator p on a Hilbert 
space Ti ~ C 2 , which can be expanded with respect to the standard Pauli basis, 
p = xo~ x + y(Ty + za z ), where I is the identity matrix and 

°* = ( I I ) , "v = ( • o* ) ' a * = ( -1 ) • (1) 

are the usual Pauli matrices, and we have * = Tr(cr ,,,/?) for * 6 {x,y,z}. Any single 
qubit gate can be written 

R A (a) = exp (i^ n • <x) , (2) 

where cr = ((T x ,ay,a z ), which can be interpreted as a rotation of the Bloch vector 
s = (x,y,z) € K 3 by an angle a around the axis n = (n x ,n y ,n z ). We further note 
that the norm of the Bloch vector ||s|| < 1, with equality if and only if p represents a 
pure state. 

The ability to implement universal quantum operations on n qubits requires a 
minimal set of elementary gates, which typically comprises several essential single 
qubit gates such as the identity I, the S and T gates and the Hadamard gate 

5" = i?z(|), T = R & (j), Had = i? fi (7r), (3) 

where n = ^(1,0,1), and a single universal two-qubit gate such as the CNOT 
gate [10 . Often the 7r-rotations about the x, y and z-axis, X = cr x , Y = a y , Z = a z , 
respectively, are also included for convenience. 

To protect quantum information from errors due to noise and decoherence several 
physical qubits are required to encode a logical qubit in a way that enables us to 
recover quantum information from corrupted qubits. Different types of encodings 
exists to protect against different types of errors, but the most common encodings 
are decoherence-free subspaces and stabilizer codes [TU]. In this paper we focus on 
the latter type of encoding, where errors can be corrected by performing (projective) 
syndrome measurements on encoded qubits and applying the correcting quantum gates 
to the corrupted qubit. The minimum number of physical qubits required to encode a 
single quantum bit of information such that we can recover from a single bit or phase 
flip error on any of the physical qubits is the five-qubit stabilizer code based on the 
encoding 

|0l) = ^Iiooooo) + |10010) + |01001) + |10100> 

+ |01010) - |11011) - |00110) - |11000) 
- |11101) - |00011) - |11110) - |01111) 
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- |10001) - |01100) - |10111) - |00101)] (4) 
\1 L ) = j[|imi) + 101101) + |10110) + |01011) 

+ |10101) - |00100) - |11001) - |00111) 

- |00010) - |11100) - |00001) - |10000) 

- |01110) - |10011) - |01000) - |11010)], (5) 

where 1 00 . . .) is the usual shorthand for the tensor product |0) (g) |0) ® Thus, a 

single qubit logic gate on encoded qubits is equivalent to a five-qubit gate represented 
by a 32 x 32 unitary matrix. There is a certain degree of freedom in the definition of 
logic gates but a logic X-gate, for instance, must obviously swap the logic states |0t) 
and |1l). Moreover, the key premise of fault-tolerant logic gate construction is that 
the application of a particular logic gate to a corrupted qubit, should not increase the 
number of errors. Thus, an encoded qubit corrupted by a single bit or phase flip error 
on a physical qubit must be mapped to a state with a single bit or phase flip error to 
enable subsequent error correction. If X n is denotes a bit flip applied to the nth qubit 
then the requirement of unitarity prevents us from directly correcting errors as we 
cannot map both |0l) and the corrupted state X n |0l) to |1l)- With this in mind, it 
is straightforward to derive suitable matrix representations for the fault-tolerant gates. 
For example, the logical Xl gate must map |0l) to |1l) and vice versa. Furthermore, 
an input state corrupted by a single bit flip error on the nth qubit, X n |0i) should be 
mapped to the output state X n i-e., the target state |1_l) with a bit flip on the 
nth qubit, or more generally 

X L = (0l| + J2 t X « I 1 *) (° L \ X " + XnX ™ I 1 *) (0 L \X m X n ] + c.c. (6) 

where c.c. denotes the complex conjugate. The same relationship between the input 
and output states should hold when the bit flip errors X n are replaced by phase flip 
errors Z n . This fully determines the logic X^ gate. We can similarly derive explicit 
representations for all other logic gates. Fig. [I] shows the structure of the resulting 
unitary operators corresponding to the target gates. 

Given a register of identical and non-interacting qubits, most fault-tolerant 
quantum logic gates can be implemented transversally, i.e., by applying the same 
single qubit operation in parallel to each qubit in the register. For instance, it 
is can easily be verified that the logical X-g&te defined by (JsJ) is simply Xl — 
X(E)X(g)X(E)X(g)X, i.e., an X gate applied to each physical qubit, However, even 
if we are given ideal non-interacting physical qubits, not all quantum logic gates can 
be implemented transversally. For instance, the T-gate above, which is required for 
universality, cannot be realized this way. Indeed, implementation of the T-gate for 
the five-qubit code above is generally complicated, requiring auxiliary qubits and 
quantum teleportation [10 . Furthermore, for real physical systems inhomogeneity 
and uncontrollable interactions may make simple transverse implementation of fault- 
tolerant gates impractical, if not impossible. 

3. Model systems and geometric control schemes 

Quite a few proposed realizations for quantum computing architectures can be formally 
thought of as linear arrays of two-level systems (qubits or pseudo-spin- ^-particles). In 
the ideal case one usually considers individual spins that are identical, individually 
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Figure 1. Real and imaginary parts of 32 X 32 unitary operators corresponding 
to the encoded logic gates. 



addressable and controllable, and assumes the couplings between all qubits are fully 
controllable. In practice, however, these requirements are difficult to meet as the 
coupling between spins is often not controllable, and the qubits may not be identical or 
selectively addressable. In this case simple transverse implementation of quantum logic 
gates is not possible even for simple gates such as the X-gate, and an optimal control 
approach seems the most promising way to overcome such obstacles to implement 
quantum logic gates with limited control. 

In the following we consider two control paradigms (a) global control and (b) 
limited local control, as shown in Fig. [2] In case (a) the physical qubits comprising 
the logic qubit are controlled entirely by a global field that interacts with all physical 
qubits simultaneously. In case (b) there may also be a fixed global field but the 
actual control operations are performed using local control gates, such as local control 
electrodes that tune individual spins in and out of resonance with a global field. 
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Figure 2. Spin chain controlled (a) by modulation of global fields and (b) using 
a fixed global field and local electrodes that tune individual qubits in and out 
of resonance with the field. B x (t), B y (t) and B z (t) indicate the x, y and z 
components of a globally applied field. Vj, refer to control parameters for local 
actuators, e.g., voltages applied to control electrodes. 
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Table 1. Euler decomposition of elementary gates in terms of rotations about 
x and y axis and total pulse length in units of 2f2/7T, where f2 is the effective 
(average) Rabi frequency of the pulse. 



3.1. Global Control 

If the qubits are identical and there are no couplings between them, then transversal 
logic gates can be easily implemented in the global control setting by simply finding 
a control pulse that implements the desired gate for a single qubit and applying it to 
all physical qubits. The elementary single qubit gates can be implemented, e.g., by 
applying a sequence of simple pulses effecting rotations about two orthogonal axes, 
using the Euler decomposition 

R&(0) = i? x (a)i? y (/3)i?x( 7 ) (7) 

for suitably chosen values of a, j3 and 7, as shown in Table [I] Rotations about the x- 
or y-axis can be performed by applying a pulse of the form f(t) = B(t) cos(ujt + (j>) 
with suitable frequency lu, pulse area / t B(t) dt and phase <j) 0. 

The problem becomes more complicated when the physical qubits are not exactly 
identical, i.e., when there is inhomogeneity resulting in the qubits having slightly 
different resonance frequencies, for instance. In this case we can switch to frequency- 
selective addressing, i.e., try to implement rotations on individual physical qubits by 
applying frequency-selective geometric pulse sequences in resonance with each qubit, 
either sequentially or concurrently. However, the pulse amplitudes (Rabi frequencies) 
have to be much smaller than the frequency detuning between the different qubits in 
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this case to avoid off-resonant excitation effects. For instance, consider the system 
H[f{t)] =H Q + f(t)H x with 

ffo = ^E^ n) . ffi = 5E«i n) , (8) 

71=1 71— 1 

where f(t) — ^B x (t) and (a^) denotes a five-fold tensor product, whose nth 

(2) 

factor is a z (a x ) and all others arc the identity I, e.g., 

can implement a logic X-gate by concurrently applying five Gaussian 7r-pulses 

/„(*) = qV^e-o^*-*^ 2 cosKi), (9) 

i.e., choosing f(t) — Y] ^ fn(t)i where ui n are the resonant frequencies 
{6,8,10,12,14} above. t n is the length of the nth pulse. If the pulses have equal 
lengths and are applied concurrently then t n — tp for all n. We can also apply the 
pulses sequentially but this significantly increases the total gate operation time as 
tp = ~^2 n t n in this case, and we will not consider this case here. 

We can quantify the gate fidelity modulo global phases as the overlap of the actual 
gate implemented U(tp) with the target gate W, 

\Tw{U{t F ))\ = ^\Tv[W^U{tp)]\, (10) 

where the factor -h, N = dim7Y, is a normalization factor included to ensure that \!F\ 
varies from to 1, with unit fidelity corresponding to a perfect gate. For the model 
system above simulations suggest that we require at least approximately tp = 440 
time units for Gaussian pulses with q — 0.01 to achieve > 99.99% fidelity. 

Fig. [3] shows that the frequency-selective pulses perform as expected if the qubits 
are non-interacting. The control performs poorly, however, if non-zero couplings 
between adjacent physical qubits are present, e.g., if we add a simple uniform Ising 
coupling term 

4 
n=l 

between neighboring qubits, then the fidelity drops from 99.99% for J = to 99.78% 
for J = 0.0001, to 86.90% for J = 0.001, to 26.97% for J = 0.01. Thus, even if the Ising 
coupling frequency J is four orders of magnitude smaller than the median transition 
frequency tu — 10, the simple geometric pulse sequence becomes ineffective. This can 
be attributed to the creation of entanglement between the physical qubits over the 
duration of the pulse. The evolution of the five physical qubits on the Bloch sphere 
under the five concurrent 7r-pulses in Fig. [4] shows that the presence of the couplings 
completely changes the dynamics. In the absence of couplings the composite Gaussian 
pulse achieves the desired simultaneous bit flip for all the qubits, but in the presence 
of even moderate Ising couplings between adjacent qubits the pulse sequence becomes 
ineffective. A plot of the length of the Bloch vectors (Fig. |5| for J = 0.01 shows a 
marked decrease for all qubits, indicating the development of significant entanglement 
between the qubits. Due to the nature of the couplings (Ising-type) the development 
of entanglement is not immediate. If we start in the product state 1 00000) then the 
effect of the couplings is felt only after sufficient coherence has been created by the 
pulse. The product state 1 00000) was chosen as an initial state, although it is not a 
logic state of the stabilizer code, because it serves well to illustrate the dynamics of 
the driving system. 
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Figure 3. A field composed of five concurrent Gaussian 7r-pulses with frequencies 
{6, 8, 10, 12, 14} achieves a logic X-gate in the absence of interqubit coupling 
(J = 0) but it is ineffective, for J = 0.01, achieving a maximum fidelity 
l^l of only 27%. The high-frequency oscillations of the fields and fidelity are 
due to both being computed and shown in the stationary laboratory frame (no 
approximations) rather than in a multiply-rotating frame. 
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Figure 4. Trajectories (red) traced out by the Bloch vectors of the individual 
qubits on the Bloch sphere subject to the five concurrent 7r-pulses above for J = 
(top) and J = 0.01 (bottom). (• initial state, -Jt final state, I target final state 
(south pole) 
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Figure 5. Length of the Bloch vectors s' n '(t) of the five physical qubits as 
a function of time for J = 0.01. Its decrease is evidence of the formation of 
entanglement between the individual qubits. 
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3.2. Local control 

The second control paradigm we consider local control through selective detuning of 
individual dots in a quantum register via local voltage gates, as recently considered, 
e.g., in [IT]. The individual qubits here are subject to a fixed, global magnetic 
field B(t) as well as local voltage gates that allows us to control the energy level 
splitting oj n of each of the qubits. If there is no coupling between qubits, setting 
B(f) = B(cosujt, — sinwi, B z ) leads to a Hamiltonian of the form 

H = X E lnB[cos{ut)o-^ - S in(„i)4™)] + <"*CV(f))<tf° (12) 

n n 

where uj n (Y) indicates that the dependence of the energy level splitting oj n of the 
nth on external control voltages V. It is convenient in this case to transform to a 
moving frame Uo(t) = exp(i£w| ^ n o~i ) rotating at the field frequency u, in which 
the Hamiltonian takes the form 

n n 

if we set ft = \B^ and u n {t) = ±Aw„(V(t)) with Aw„(V(t)) = w„(V(i)) - w, 
and for simplicity assume 7„ = 70 for all n. In the following we do not consider 
the architecture-specific functional dependence of resonant frequencies on the gate 
voltages V and cross-talk issues [T3J [T5] and take u(t) = (ui(t) , u 2 (t) , . . . , u 5 (t)) to be 
independent controls. 

In the absence of interactions, it is again relatively straightforward to implement 
arbitrary rotations on any qubit using a simple geometric control design. We can 
clearly perform simultaneous rotations about the x-axis on all physical qubits simply 
by choosing the detuning parameters u n (t) — 0. To rotate qubits about the y-axis we 
can make use of the following result [11 

Ry{A<P) = R^)R^)R^)R^), (14) 

where R^(tt) is a tt rotation about the x-axis and R&(tt) is a n rotation about 

n = cos </>x + sin <fii, (15) 

where cos0 = Q/Q and sin0 = u/Qq with fi = yf2 2 + u 2 and u n (t) = u. If we can 
achieve detunings of the same magnitude as the fixed coupling, u = f2, then 0=5 
and we can therefore implement 7r-rotations about the axis n = (1, 0, 1) T , which 
corresponds to a Hadamard gate, and rotations about the z-axis using 

R ± (6) = Hadi? i (6»)Had. (16) 

Table |2] summarizes the sequences of rotations that implement the elementary single 
qubit gates, and Fig. [6] shows the evolution of the components of the Bloch vector 
for each of gates assuming piecewise constant controls. The transversal gates can 
be implemented by simply applying the same sequence of rotations on each qubit in 
parallel. 

An advantage of the local control scheme compared to frequency selective global 
control is potentially much faster gate operations because the gate operation times are 
determined by the Rabi frequency O of the global field, which is not limited by the 
need to avoid off-resonant excitation, although the need to be able to induce detunings 
on the order of imposes some constraint in practice. An added bonus of the shorter 
gate operations times is that the scheme is far less sensitive to the presence of small 
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Figure 6. Time-evolution of Bloch vector components s = (a;, z) for various 
elementary gates using local geometric control with time in units of where 
Q is the Rabi frequency of the fixed global field. 



couplings between the qubits. Adding a fixed Ising coupling term of the form (111 
does reduce the gate fidelities, but not nearly as much as was the case for the global 
control scheme as the much shorter gate operation times mean that the qubits have 
far less time to interact with each other. 

A significant disadvantage of this control scheme, however, is that the 
implementation times for different gates using the type of simple geometric pulse 
sequences described differ significantly, which is problematic in a scalable architecture 
where local operations on different logical units are to be implemented concurrently. 
Although the implementation times in units of the Rabi frequency for different gates 
also varied in the global control case, in the former case the Rabi frequencies of 
the pulses were variable (within a certain range), allowing adjustment of the gate 
implementation times. This adjustment is not easily possible in the local control case 
because the Rabi frequency f2 of the global field must remain fixed. 

In summary, we can in principle implement quantum logic gates on encoded qubits 
using simple control pulse sequences based on geometric control ideas in both the local 
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and global control setting. However, the simple schemes have serious shortcomings. 
Frequency-selective global control schemes, for instance, requires comparatively long 
pulses to minimize off-resonant excitation. In the local control setting it is difficult 
to implement different gates concurrently on different qubits due to variable gate 
operation times. Both schemes also assume non-interacting physical qubits, and tend 
to perform poorly in the presence of fixed couplings due the creation of entanglement 
between qubits. Also, some essential gates such as the T-gate cannot be implemented 
transversally at all. This raises the question if we can overcome such problems using 
optimal control, which we will consider in the next sections. 



4. Optimal Control 

The problem of finding the control pulses that produce the desired quantum gates can 
be treated as an optimization problem. The procedure in general involves choosing 
a measure of the gate fidelity to be maximized, suitably parameterizing the control 
fields and finding a solution to the resulting optimization problem. If the evolution of 
the system is given by a unitary operator U(t) satisfying the Schroedinger equation 



M 

H o+^2 u m (t)H„ 

in— 1 



U(t,t ), (17) 



where Hq is a drift Hamiltonian and H m , m — 1, . . . M < oo, are control Hamiltonians, 
the task is to find a control u(t) that maximizes the gate fidelity 

Tw[U{t F M)] = jjZtTc[W*U{t F ,t )] (18) 

or similar figure of merit for some fixed target time tp [J] 

To derive an iterative procedure for finding optimal controls, note that changing 
a given control u by some amount Au changes the corresponding propagator (see 



Appendix A ) 



A[/ u (Mo) = tWu(Mo) - U u (t,to) 



'to 

The corresponding change in the fidelity is 

AT{t F ) = T u+ Au(t F ) - FJJLf) 



[ ^Au m (T)lf u (£,T)ff m J7 u+ Au(T,i )dT. (19) 
Jtn m 



^Tr[W^U u+ ^ u (t F ,t ) - U u (t F ,t ))] (20) 
2 

Nh- 



[ Au m (T)^Tr[W^U u (t F ,T)H m U u+Au (T,t )}dT. 



| The definition of the fidelity \l8\ is stricter than the definition of ^-^[11 (tp, trj)]| in \10\ , which 
disregards the global phase of the target operator. I n g eneral, one does not care about the global 
phase of the operator U, whence the weaker measure | |10| l we used earlier is sufficient. Although, one 
can modify the following derivation to obtain an explicit algorithms to maximize the [^(^({7)[, the 
expressions are more complicated due to the absolute value, which is why we focus on optimizing the 
stricter fidelity measure FyylU), which takes the global phase into account here. One could argue 
that this is slightly unfair to the optimal control algorithm as we have actually made the problem 
harder. However, since we are working in 5(7(32) the only phase factors we are excluding are roots 
of unity. Although these can be problematic, they did not seem to pose a problem for the algorithm 
in our case, i.e., we had no problem finding solutions with W^U close to the identity I, which is why 
we did not consider modifying the algorithm to optimize ^^((J)! instead, although this could be 
done if necessary. 
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This shows that setting 

Au m (t) = e m (t)%Tr[W^U u (t F ,t)H m U u+ Au(t,t )} (21) 

with e m (t) > for all t will increase the fidelity at the target time t F . Thus, a basic 
algorithm for maximizing T is to start with some initial guess u (i) and solve the 
Schroedinger equation iteratively while updating the control according to the rule 

u'" +1 »=u<"'+Au'"), (22) 



with Au'"' as in (21). Note that this is an implicit update rule as the RHS of (21 1 



depends on Au(t), and the search space is infinite dimensional, consisting of functions 
on defined on an interval of the real line. 

To derive useful practical algorithms the controls u m (t) need to be discretized, 
and a simplified explicit update rule is desirable. The simplest and most common 
approach to discretization is to approximate the continuous fields by piecewise 
constant functions, i.e., we divide the total time t F into K steps, usually, though 
not necessarily, of equal duration At — t F /K, and take the control amplitudes to be 
constant during each interval Ik — [tk-i,tk)- We then have 



Fw\U(t F )] = -5JTr 



wiuWuF-v-'-uPuP 



where the step propagator for the fcth step is given by 

M 



exp 



•rAt(flo+5Z 

" /;//»• 



(23) 



(24) 



and Umk is the amplitude of the mth control field during the fcth step. Furthermore, 
if we change the control field by Au m {t) — Au m k for t G Ik then the change in the 



fidelity at the target time t F is (see Appendix B) 
2 



AT{t F ) 



Au mk ^s Tr 



Jj{k+l) \Jj(k)jj{k-l) 



. u. 



(1) 

u+Au 



(25) 



with Uu and (7^ +4u as defined in (24 1 and 

rt k 

AC/W = / U u (t k , T)H m U u+ *u(T, t fc _i)] dr. 

This shows that the fidelity will increase at time t F if we change the field amplitude 
u m k in the time interval Ik by 



(26) 



Au mk = Cmfe^Tr 



W^U^---U^AU^(t)U^l...U^ 



Au 



for e m k > 0. This is still an implicit update rule as AU^ on the RHS depends 
Au rn k- However, when At = tk — ifc-i and Au m k are sufficiently small then we can 



approximate AU, 
update rule 



O) 

m 7 



e.g., setting AU 



(fc) 



U u m. -i 



which yields the familiar explicit 



Au ri 



■ U 



(1) 

u+Au 



(27) 



In practice we can now solve the optimization problem by starting with an initial guess 
for the controls ufj' , calculating U^\t F ) and then integrating alternatingly backwards 
and forwards while updating the controls in each time step according to the explicit 



update rule (27 1 until no further improvement in the fidelity is possible. Ignoring 



errors to due numerics (e.g. limited precision floating point arithmetic, etc), if the 
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algorithm is implemented strictly as presented and the e m k are chosen in each time 
step and iteration such that the value of the fidelity at the final time does not decrease, 
then convergence is guaranteed as the fidelity is nondecreasing and bounded above. 
However, we cannot guarantee that the fidelity will converge to the global maximum, 
especially when there are constraints on the field amplitudes or time resolution At of 
the fields. The algorithm will stop when it is no longer possible to change u m k for any 
m or k such as to increase the final fidelity. 



The update rule (27 1 is similar to the update rule of the familiar GRAPE 
algorithm [Tl], but in the GRAPE algorithm the update is global, i.e., the field 

amplitudes u m k for all times tk are updated concurrently in each iteration step based 

(k) 

on the propagators J7u using the fields from the previous iteration step, which is 
equivalent to setting 



Au m , = e m 3 Tr W*U™ ■ ■ ■ U^H m U^ . . . 



(28) 



for all m and k. This global update has certain advantages such as easy 
parallelizability, but it requires the field changes in each step to be small to maintain 
the validity of the approximations involved and ensure an overall increase in the 
fidelity The parameter e m is critical in this regard and must be chosen carefully to 
ensure monotonic convergence. One way this can be achieved is by choosing e m to be 
very small (steepest decent) but this comes at the expense of very slow convergence. In 
practice one would therefore usually employ step-size control in the form of lincscarch, 
quasi-Newton methods or conjugate gradients to accelerate convergence. 

Despite the similarities in the update rules, the sequential local update method 
is quite different. With sequential local update the fidelity increases in every time 
step, not just for every iteration, and we can make much larger changes to the field(s) 
in each step, with the magnitude of the allowed changes limited mainly by explicit 
constraints on the field amplitudes and the approximations involved in transforming 
the implicit update rule to an explicit one, which could be avoided or improved, and 
discretization errors. The sequential local update rule also allows us to easily enforce 
certain local constraints. For example, a constraint on the field amplitudes of the form 
|w m (f)| < C can be trivially incorporated as we simply have to choose t m k m each 
time step k so that \u m k + Au m k\ < C. Assuming the initial controls is choosen such 
that the constraints are satisfied, \u m k\ < C, this is always possible, if necessary by 
setting e mn — for this time step. In this case there will, of course, be no increase 
of the final fidelity during this time step but we will generally be able to continue to 
increase the fidelity in the next time steps, although such an amplitude constraint may 
eventually prevent us from increasing the fidelity further for all time steps, in which 
case we could be left with a solution for which final fidelity is less than its optimum 
value. Of course, amplitude constraints can also be incorporated in global update 
schemes, albeit not quite as trivially. Another perhaps more interesting feature of the 
sequential update method is that it enables us use variable time steps to increase the 
fidelity further should a given time resolution of the fields prove insufficient to achieve 
high-enough fidelities. 

The sequential local update method is similar to a class of methods often referred 
to as the Krotov method, based on work by Krotov [T5J [16], adapted to quantum 
systems by Tannor et al |17j , although our formulation is for unitary operators rather 
than quantum states and there is no penalty on field energy. 
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5. Optimal Control Implementation of Encoded Logic Gates 

We now apply optimal control algorithm described in the previous section to the 
problem of implementing fault-tolerant logic gates for the five-qubit stabilizer code. 
In particular, we wish to implement a full set of elementary logic gates in a fixed 
amount of time for systems subject to both inhomogcncity and interactions between 
adjacent qubits. For the global control system we choose the total Hamiltonian 

ffte)[u(i)] = ^ + •'E 4"M" +1) + «i(*) E 4° + u ^ E 4 n) ( 29 ) 

n— 1 n—1 n—1 n—1 

with uj n £ {6, 8, 10, 12, 14}, and the aim is to optimize the two fields ui(t) and v,2(t), 
which roughly correspond to the x and y components of an external electromagnetic 
field, to implement the desired encoded logic gates. For the local control case we 
choose the Hamiltonian 

[ U (t)] = n e 4*° + J E 4 n) 4 n+1) + E u ™ (*)4 W) (3°) 

71 n—1 n 

where u n (t) for n = 1,2,3,4,5 are controls that induces local detunings and fi = 10 
is the Rabi frequency of the fixed external coupling field. Here we have chosen units 
such that H = 1 . 

In order to find the optimal solutions using the algorithm outlined above we must 
choose the number of time-steps K and the total time tp and a suitable initial control 
u(°). To ensure that the gates are fast and do not require controls with extremely 
high time resolution, we aspire to find solutions for tp and K small. Unfortunately, 
there are no simple rules for choosing tp and K. With some experimentation we were 
able to find controls that achieved gate fidelities of approximately 99.99 % for the 
global control example with J = 1 for tp = 125 with K = 1250 time steps. For the 
local control example tp = 30 and K — 300 proved sufficient. The gate operation 
times achieved using optimal controls in the global control setting, although much 
longer than for the local controls, are still favorable compared to frequency-selective 
geometric control pulses, for which even the A-gate required approximately 500 time 
units for J = and the method failed for J greater than 10~ 4 . Although the local 
gates for 17 = 10 are much slower than what is theoretically possible using local 
geometric control in the absence of coupling, the optimal controls have the distinct 
advantage that they enable us to implement all gates in a fixed amount of time and 
deal with non-zero J-couplings. Also, the implementation of the T-gate, which can 
not be implemented transversally, presents no greater challenge for optimal control 
than any of the other gates. In fact, the simulations suggest that the most difficult 
gate to implement in the global control setting is the legate. In the local control 
setting there appeared to be little difference in the difficulty of implementing different 
gates. 

For comparison with the geometric control results, Fig. [7] shows the trajectories 
of the physical qubits subject to the optimal controls for the logic A-gate. Again, 
the qubits were assumed to be initialized in the product state 1 00000). This is not 
a logic state but allows us to easily check visually that the gate effects the desired 
simultaneous bit flip on all physical qubits even in the presence of significant intcr- 
qubit couplings (J — 1). The single-qubit trajectories start out at the surface (north 
pole) of the Bloch sphere, corresponding to the product state 1 00000). They descend 
inside the Bloch ball as the qubits then become entangled, and the individual qubits 
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Qubitl Qubit2 Qubit 3 Qubit 4 Qubit 5 




Figure 7. Evolution of physical qubits initialized in the state 1 00000} subject to 
the optimal controls required to implement the X-gate. 
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Figure 8. Controlled evolution of the initial state |*(0)> = ([Ox,) + \1 L ))/V2 
projected onto the 2D subspace spanned by |0x) and |lx) for different logic gates 
for optimal control examples with global and local control. 



follow different paths, but at the final time the trajectories of all physical qubits 
converge to the south pole of the Bloch sphere, as desired. 

To get an idea of how the optimal controls perform for the other gates beyond 
the gate fidelity, Fig. [8] shows the projection of the controlled evolution onto the two- 
dimensional subspace of the Hilbert space spanned by the logic states \0l) and 
for the initial state |^(0)) = ^(\0l) + The Bloch vector was defined to be 

s x {t)=2$t[(0 L | *(*))<*(*) | 1 L >] (31) 

S „(f) = -29f[(Qi | | li)] (32) 

Sz (t) = \(l L \*(t))\ 2 -\(0 L \*(t))\ 2 , (33) 

where |W(£)) = U(t) |^(0)) and U(t) is the full 32 x 32 unitary operator describing 
the evolution of the system under the given optimal control. The plots show that the 
evolution of the system is very complicated. Although the initial state is in the 2D 
subspace spanned by the logic states |0l) and |1l), the fact that the trajectories (red) 
are concentrated in the interior of the Bloch ball indicates that the system spends most 
of the time outside the two-dimensional logic subspace, but returns (mostly) to this 
subspace near the final time. Although the gate fidelities as defined in (18 1 are about 
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Figure 9. Example of optimal control solution to implement the non-transversal 
T-gate for the five qubit stabilizer code using two independent global fields for 
system (|29|. 



99.99 %, there are small but noticeable differences between the actual final states (-fc) 
and the final states (blue squares) for an ideal gate. This is due to the fact that the 
normalized gate fidelity for high-dimensional gates is not a very strict measure as 



\\w- u(t F )h = - u)Hw - u)} 

= y/2N[l - F W (U)], 



(34) 



with Tw(U) as defined in Eq. (18), i.e., the distance between the operators with 



respect to the Hilbert-Schmidt norm is significantly larger than the gate error defined 
in terms of the fidelity 1 - T w [U{t F )]. For N = 32 a gate fidelity of > 99.99 % allows 
a Hilbert-Schmidt distance between the operators up to 0.08. The distance between 
the operators with respect to the operator norm ||X|| (given by the largest singular 
value of X) is usually smaller than the Hilbert Schmidt distance but still significantly 
larger than the fidelity error. A summary of different gate error measures for all gates 
is given in Table [3j 

Examples of the actual control fields obtained in both the local and global control 
setting are shown in Figures [9] and 10 for the T-gate. The control fields appear 
complicated, especially in the global control case. The local control fields appear to 
be better behaved but this is is in part due to the scaling of the time axis — the global 
control pulses are much longer and thus appear much more compressed when scaled to 
fit the width of the text. Nonetheless the controls are complicated and have a relatively 
high bandwidth. To some extent this is unsurprising considering the complexity of the 
system, control goals, and constraints imposed on the controls in terms of pulse lengths 
and time resolutions. However, there are many possible controls that result in the same 
gate, corresponding to alternative trajectories in SU(N) connecting the identity and 
the target gate. A limited number of numerical experiments performed with different 
choices for the initial fields suggest that these lead to different solutions for the controls, 
but generally the solutions appear to have similar complexities. However, it may be 
possible to exploit the non-uniqueness of the solutions to find better solutions, cither 
by changing the parameterization of the fields or imposing carefully chosen additional 
penalty terms, which we will consider in future work. 
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Figure 10. Example of optimal control solution to implement the non-transversal 
T-gate for the five qubit stabilizer code using five independent local detuning fields 
for system i30\. 
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Had 


global 




0.9999 


0.9999 


0.9996 


0.9999 


0.9999 


0.9999 


0.9999 




\\w-u\\ 


0.0267 


0.0268 


0.0573 


0.0259 


0.0255 


0.0261 


0.0266 




\\w-u\\ 2 


0.0800 


0.0800 


0.1644 


0.0800 


0.0800 


0.0800 


0.0800 




\\W-U\\ max 


0.0061 


0.0069 


0.0163 


0.0074 


0.0064 


0.0071 


0.0065 


local 




0.9999 


0.9999 


0.9999 


0.9999 


0.9999 


0.9999 


0.9999 




\\w-u\\ 


0.0268 


0.0266 


0.0266 


0.0253 


0.0252 


0.0250 


0.0263 




\\W-U\\ 2 


0.0799 


0.0798 


0.0799 


0.0798 


0.0798 


0.0798 


0.0799 




||W-£/|| max 


0.0056 


0.0069 


0.0069 


0.0058 


0.0069 


0.0068 


0.0066 



Table 3. Different performance measures for the gates implemented using optimal 
control: Fidelity, operator distance \\W — U\\ (largest singular value of W — U), 
Hilbert-Schmidt distance \\W - U\\ 2 = ^Tr[(VK - U)i(W - U)], and maximum 
distance between matrix elements \\W — (7|| ma x = max mi „ |W ran — U m n\, where 
U = U(tp) and W is the target gate. 



The presence of fixed couplings between the qubits, which was highly detrimental 
for the geometric control schemes, does not appear to be a problem in the optimal 
control case, and indeed some of our simulations suggested that larger J-couplings 
may in fact make it easier to find optimal control solutions. This apparently counter- 
intuitive behavior prompted us to consider the effect of varying the interqubit coupling 
strength J and the detunings Aw. We can get an estimate of the difficulty of finding 
solutions by setting tp and K to constant values and comparing the number of 
iterations required to archive 99.99% fidelity for various values of J and Au>. A 
small number of iterations generally implies that it is easy to find a solution and that 
solutions may exist for significantly shorter times tF- To reduce the computational 
overhead, we considered the same model systems but with only three instead of five 
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Figure 11. Number of iterations required to find solutions with Tw{U) > 0.9999 
for (a) t F = 10, M = 80 and Aw = 2; (b) t F = 10, M = 80 and J = 2. 



qubits, and taking as target gates the standard gates for the three-qubit bit-flip code 
with code words = |000) and — |H1). For this system the run time of the 
optimization routine reduces from several hours for the five qubit code to a few minutes 
on a single 3GHz CPU. It proved easy to find solutions for the system subject to global 
control with J = 1 for tp — 10 and K = 80, so we selected these values as defaults. 



Fig. lUa) shows the number of iterations for different values of J. Although the 



number of iterations required for convergence maybe an oversimplified measure of the 
complexity, the data indeed suggests that increasing J makes it easier for the algorithm 
to find solutions, suggesting that optimal control allows us to exploit couplings that 
are usually considered undesirable to speed up gate operations, especially in the global 
control setting. The effect of changing Auj is more subtle. If Aui is very small then 
the selectivity of the pulses is very weak and the algorithms struggles to find solutions 
unless tp is increased. If Aw is very large then it is also hard to find the solutions, 
partly because the system can not use off-resonant excitations as much, but more 
likely because the bandwidth of the fields available may not be sufficient to cover the 
required frequency range. 



6. Conclusion 



We have considered the problem of implementing quantum logic gates acting on 
encoded qubits, As a particular example, the five-qubit stabilizer code was considered. 
Despite the complicated encoding, most elementary logic gates can in principle be 
implemented fault-tolerantly by applying the same local gates to all physical qubits. 
Thus, it is usually argued that the implementation of logic gates on encoded qubits can 
be reduced to the implementation of single and two-qubit gates on physical qubits, 
and that the latter can in principle be implemented using simple pulse sequences 
based geometric control ideas. However, there are several problems with this. Not 
all gates can be implemented transversally. For the five-qubit code, for instance, it is 
not possible to implement a T-gate this way, which is an essential gate for universal 
quantum computation. Moreover, for realistic physical systems inhomogeneity and 
often uncontrollable interactions between nearby physical qubits pose significant 
problems, rendering simple gate implementation schemes based transversal application 
of single-qubit geometric pulse sequences ineffective. 

Using an optimal control approach we can in principle overcome these problems 
and find effective control pulse sequences to implement encoded logic gates even for 
imperfect systems subject to limited control, e.g., if we are unable to selectively address 
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individual qubits, or switch off unwanted interactions between qubits. Perhaps even 
more importantly, the implementation of gates that are difficult to implement using 
conventional techniques such as the T-gate for the five-qubit stabilizer code, appears 
to present no greater challenge for optimal control than the implementation of any 
other gate. A potential drawback at present is the complexity of the optimal control 
solutions. However, the optimal control solutions are not unique, and the solutions 
found here should therefore be regarded more as a demonstration of principle, i.e., that 
is is possible to find solutions to the optimal control problem. Simpler solutions are 
likely to exist at least for certain problems, and future work exploring whether such 
solutions can be found by considering different field parameterizations or imposing 
suitable penalty terms that steer the algorithm away from complicated fields towards 
simpler solutions (if they exist) could be fruitful. 
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Appendix A. Derivation of Eq. (19) 

Letting for convenience 

H' {t)=H + J2 U rn(t)H m 

m 

H[(t) = J2 Au m(t)H m . 

m 

we can rewrite the Schrodinger equation for £/ u (i, to) and £/u+Au(£j£o) as 
U u (t,t ) = - ~flS(i)0u(t,to), 
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t>u+Au(Mo) = - j|[#o(*) + ^(*Pu+Au(Mo). 

Inserting an interaction picture decomposition 

t/ u +Au(Mo) = [/u(Mo)CMMo) 
into the Schrodinger equation we obtain using the product rule 

Ui(t,t Q ) = -J[[/ u (M )ti^(i)[/ u (Mo)]^/(Mo). 

n 

Rewriting this differential equation in integral form gives 



(A.l) 



U T (t,t )=l-- / [U n (T,1b)iH[(T)U n (T,h)]Ul(T,to)dT 
J to 

Multiplying both sides by U u (t 7 t ) and noting that U u (t,t ) does not depend on r 
and thus can be moved inside the integral now yields 

C/ u +au(Mo) =U u (t,t )-- / UjJt^UJj^yH^U^r^Ui^t^dT 



= U n {t,t )-- \ U u (t,T)H[(T)U u+ ^ u (r,t a )dT 



U u (t,t ) 



/ Au m (t)U u (t, T)H m U u+ A. u (T, t ) dr 

' "to m 



where we have used U u (t, to)U u (T, t o y = U u (t,T) we as well as Eq. (A.l I and the 
definition of H[ (t) . 



Appendix B. Derivation of Eq. (25) 

If Au m (t) — except for t G Ik = [tk-i,tk] then we have 
AT(t F ) = j^-Yl I F ^u m (T)^Tr[W^U u (t F ,T)H m U u+Au (T,t )]dT 



Nh 
2 

Nh- 



Y, f " Au m (T)%Tr[W^U u (t Fl T)H m U u+ ± u (T,t Q )}dT 



Furthermore, we have U u (t F ,r) = U u {t F ,tk)U u (tk,T) as well as t/ u +Au(i", ^o) = 
E/ u +Au(V, ife-i)C/ u +Au(ife-i,io) and because U u (t F ,t k ) and f7 u+ Au(ife-i,*o) are 
outside of the domain of integration we can re- write the expression for AT(t F ) as 

AT{t F ) = -|-^A Mm ( T ) 5Tr \w^U u (t F ,tk)AUWu u+ ± u (tk,t ) 



with Auln ] as in (|26[) . Finally noting that 



U u (t F ,t k ) =U^---U^ 

fc-1) TT (1) 

+Au ■ ■ ■ ^u+Au 



U u+ & u (t k , t ) — u^ + ^ u . . . 



with and U^2/\ u as defined in (24 1 gives 



